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Abstract. We propose a method for controllable preparation and detection of 
interaction-induced ferromagnetism in ultracold fermionic atoms loaded in optical 
superlattices. First, we discuss how to probe and control Nagaoka ferromagnetism in 
an array of isolated plaquettes (four lattice sites arranged in a square). Next, wc allow 
for weak interplaquette tunneling. Since ferromagnetism is unstable in the presence 
of weak interplaquette couplings, we propose to mediate long-range ferromagnetic 
correlations via double-exchange processes by exciting atoms to the first vibrational 
band. We calculate the phase diagram of the two-band plaquette array and discuss 
conditions for the stability and robustness of the ferromagnetic phases in this system. 
Experimental implementations of the proposed schemes are discussed. 
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1. Introduction 

The origin of ferromagnetism in itinerant electron systems remains an important open 
problem in condensed matter physics. Mean field approaches, such as the Hartree- 
Fock approximation [T] and the Stoner criterion [2] for ferromagnetic instabilities, 
are extremely unreliable since they overestimate the stability of the magnetic-ordered 
phases [3]. The only rigorous example of ferromagnetism in the generic Hubbard 
model [Ij, predicted by Nagaoka in 1965 [5J, was proven for a system with one fewer 
electron than half filling (i.e., one hole) in the limit of infinite interactions. Even though 
such a ferromagnetic state is an iconic example of a strongly correlated many-body state, 
it is highly unstable and counter examples indicating the absence of ferromagnetism with 
two or more holes have been found [0]. 

The experimental observation of Nagaoka ferromagnetism (NF) is a challenging 
task, as it requires a system with a finite and controllable number of holes. Even 
though there have been recent attempts to explore Nagaoka ferromagnetism using 
arrays of quantum dots [7j, the exponential sensitivity of the tunneling rates to the 
interdot distance and the random magnetic field fluctuations induced by the nuclear 
spin background have prevented its direct experimental observation. 

Here, we propose to use cold fermionic atoms in optical superlattices for the 
controllable observation of interaction- induced ferromagnetism. First, we show how 
to probe the onset of NF in an array of isolated plaquettes (four lattice sites in a square 
geometry). Next, we discuss how to engineer long-range ferromagnetic correlations, 
starting from the isolated plaquette arrays. Since weak coupling of the plaquettes 
destroys ferromagnetism [S], we instead propose to use additional atoms loaded in 
excited bands. The underlying idea is to use Hund's rule couplings [9J to favor local 
ferromagnetic alignment among the atoms in different bands and a double-exchange 
mechanism (tunneling-induced alignment of the spins) [lO] to stabilize ferromagnetic 
correlations between adjacent plaquettes. Exact numerical calculations for an array of 
weakly coupled plaquettes confirm the existence of stable ferromagnetic order in this 
two-band setup. Finally, we discuss methods for experimental preparation and detection 
of the ferromagnetic correlations. 

2. Probing Nagaoka ferromagnetism in superlattices. 

The low-energy physics of fermionic atoms loaded in the lowest vibrational band of an 
optical lattice is well described by the Hubbard Hamiltonian: 



where Jry = J is the tunneling energy, and U is the onsite interaction energy. In 
Eq. ([l]), Cj.^ are fermionic annihilation operators, = c\^c^^ are number operators, 
r = 1, . . . L labels the lattice sites, and (r, r') in the summation indicates that the sum 
is restricted to nearest neighbors. 




(1) 



(r,r'),CT 
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Before starting let us explicitly state Nagaoka theorem [5] : "Let the tunneling matrix 
element between lattice sites r and r' be negative, J„i < 0, for any r ^ r' and U = oo 
and let the number of fermions be N = L — 1, with L the total number of sites. If 
the lattice satisfies certain connectivity condition, then the ground state has total spin 
S = N/2 and it is unique, apart from a trivial {N + l)-fold spin degeneracy". 

The notion of "connectivity" requires that each site in the lattice is contained in 
a loop (of non- vanishing J^r') and furthermore that the loops should pass through no 
more than four sites [H]. The requirement Jrr' < is just the opposite of that assumed 
in most discussions of the Hubbard model. However, in bipartite lattices, there is always 
a canonical transformation connecting J^r' < and Jrr' > 0. 

From the conditions of the theorem it is clear that the minimal geometry to observe 
Nagaoka ferromagnetism is a triangle, however a triangle is a trivial example since in 
this case either the ground state is always a singlet ( case Jr,r' = J > ) or it is always 
a triplet ( J < 0) [|j The first non trivial example of Nagaoka crossing takes place in a 
plaquette loaded with three fermions (Fig. 1(a)). 




Figure 1. (a) Schematic representation of a plaquette. (b) Low energy spectrum of 
three fermions. Sohd hnes and dash-dotted are 5 = 1/2 while dashed curve corresponds 
to S — 3/2. Inset: Zoom in of the spectrum at the Nagaoka crossing. Solid lines are the 
spectrum for zero gradient field while dashed lines are the spectrum at finite gradient 
field. 

The energy levels of a plaquette loaded with three fermions can be classified 
according to the total spin S and the symmetries of the wave function. It is known (e.g., 
Ref. [8]) that for U < Ut ^ 18. 6J the ground state is a degenerate doublet 5* = 1/2 
state with r = ± ipy symmetry (the wave function changes phase by ±7r/2 upon 7r/2 
rotation). For U > Ut ^ 18.6 J, the ground state becomes a ferromagmetic S = 3/2 
state, in agreement with the Nagaoka theorem (Fig. [l]^b)). We denote these eigenstates 
as 15* = 1/2, Sz,T = ±) and \S = 3/2, S^) with = —S,...S and recall that the 
energies are independent of the 5*2 value. The onset of Nagaoka ferromagnetism can be 
understood as competition between the kinetic energy and superexchange interactions. 
In the f/ — )■ oo limit, double occupancies are energetically suppressed, and the low-energy 



t Since a triangular lattice is not a bipartite lattice, Nagaoka theorem only holds for J < 0. 
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states are singly occupied with an energy spectrum given by E = ±2 J, ±73 J, ± J, 0. The 
relevant low-lying eigenstates are the ones with E^^^/"^ = —2J and £'•^=1/2 = —y/3J. 
As U become finite, while the fully polarized states remain eigenstates for any U and 
their energy is unaffected by interactions, the \S = 1/2,S'^,±) states acquire some 
admixture of double occupancies, which tend to lower their energy. The energy shift in 
the 5 = 1/2 states can be calculated by using second order perturbation theory, yielding 
^5=1/2 _ _y/3j _ 5j_^ rjj-^g Nagaoka crossing occurs at the Ut/J value when the two 
energies become equal, Ut = 5/(2 — \/3)J ~ 18. 66 J, in very good agreement with the 
exact diagonalization. 

An array of plaquettes can be created by superimposing two orthogonal optical 
superlattices formed by two independent sinusoidal potentials that differ in periodicity 
by a factor of two, i.e., V{x) = K/2 cos(47rx/As) — VJ/S cos(27rx/As), where Vi is the 
long lattice depth, Vg is the short lattice depth, and Xg is the short lattice wavelength. 
By controlling the lattice intensities, it is possible to tune the intra- and interplaquette 
tunneling and, in particular, to make the plaquettes independent. Here the axial optical 
lattice is assumed to be deep enough to freeze any axial dynamics. To load the plaquettes 
with three atoms, one can start by preparing a Mott insulator with filling factor three in 
a 3D lattice and then slowly split the wells along x and y. Since only two fermions with 
opposite spin can occupy the lowest vibrational level, loading three fermions per site 
requires populating the first excited vibrational state before splitting the wells. This 
loading procedure creates plaquettes with S = 1/2. 

Since S and Sz are conserved quantum numbers in clean cold atom set-ups, to 
probe the Nagaoka transition we require the presence of a weak magnetic- field gradient. 
We choose for this case a field pointing along z with a constant gradient along the x 
direction B(x) = ^^|^z, where /ib is the Bohr magneton g is the gyromagnetic factor. 
The magnetic-field gradient couples the |3/2) state with some linear combination of 
|l/2, l/2,±) which we denote as |l/2, 1/2,1), through a Hamiltonian matrix element 
-^3/2,1/2 = —2/3(1 + \/3)SEb and leaves another linear combination of |l/2, 1/2, ±), 
which we denote as 1 1/2, 1/2,2), uncoupled. It consequently transforms the crossing at 
Ut into an avoided crossing [see inset of Fig. [l] (b)] which can be used to adiabatically 
transform the |l/2, 1/2, 1), ground state for U < Ut, into |3/2) as U is slowly increased. 
U can be tuned by a Feshbach resonance in the presence of an external magnetic field. 

The variation of (S) in the ground state below and above the Nagaoka point can 
be inferred from a band mapping analysis [12] after collapsing the plaquettes into single 
wells: The states with 5 = 3/2 and 5=1/2 will occupy three and two bands respectively 
when a plaquette is transformed into a single site. Since the state |l/2, 1/2,2) remains 
decoupled, one can improve the band mapping signal by adiabatically turning on the 
magnetic field gradient at the same time the plaquettes are prepared. In this way the 
initial plaquette configuration will be almost in a pure |l/2, 1/2, 1) state and we will be 
able to fully transform the state into the |3/2) at U > Ut. 

The energy difference between the |3/2) and 1 1/2, 1/2, ±) states can also be probed 
dynamically. After preparing the state \ip{0)) = cosa|3/2, 1/2) + sina|l/2, 1/2, 1) 



Probing Interaction- Induced Ferromagnetism in Optical Superlattices 



5 



one can suddenly turn off the magnetic-field gradient. By measuring the Neel order 
parameter or spin imbalance along the x direction [Ns{t) = l/2(X]r=i 2 '^t?- ''^ir ~ 
Sr=3 4 ^t'r~'"'4r)] as a fuuctiou of time, one can track the Nagaoka point by the oscillation 
period of {Ns{t)) oc cos 

[(^5=3/2 _ E^=^/'^)t/h]. As U/J approaches Ut/J, the period 
will become very long, indicating that the character of ground state has changed. This 
simple treatment ignores the admixture of states with double occupied sites in the 
|l/2, 1/2,1) state. When included, the excitations introduce fast but small oscillations 
of frequency J. Comparisons between the exact and analytic solutions are shown in 
Fig. g The spin imbalance Ns{t) can be experimentally probed by first splitting the 
plaquettes into two double wells and then following the same experimental methods used 
for measuring superexchange interactions [13] that rely on band-mapping techniques and 
a Stern-Gerlach filtering. 
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Figure 2. Normalized spin population imbalance. At the Nagaoka crossing, the 
envelope frequency becomes very long, indicating zero-energy splitting between the 
1 3/2) and |l/2) levels. 



The constant magnetic field needed for tuning a Feshbach resonance does not affect 
the dynamics since the relative energy spacing of the various levels within a plaquette is 
insensitive to such magnetic fields. The big advantage of this probing method is that it 
does not require fixing the same magnetization for the various plaquettes. Consequently, 
we can relax the temperature constraint for preparing the Mott insulator used for the 
initial loading. The insensitivity of this probing method to the initial magnetization can 
be understood by the fact that the dynamic taking place in a plaquette initially loaded 
with Sz = —1/2 is identical to that described for the 5*^ = 1/2 case. Furthermore, 
preparation of plaquettes with Sz = ±3/2 is energetically suppressed due to the large 
vibrational energy spacing. It should be stressed that finite temperature effects only 
determines the fidelity of the initial preparation of a Mott insulator. Once the Mott 
insulator is achieved, the observation of ferromagnetism depends on the efficiency of 
implementing the adiabatic manipulation discussed above. 
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3. Engineering long-range ferromagnetic correlations. 

We now study the more general case in which one allows a weak interplaquette 
tunneling, J', by lowering the long lattice depth along both the x and y directions 
(or along only x). This procedure generates a 2D (ID) array of plaquettes. In the 
Nagaoka regime {U/J > 18.6) to zero order in J', the many-body ground state has 
a degeneracy of 4^ (A^ is the number of plaquettes) and is spanned by states of the 
form \^)s^i,...s^N — Ylil^ = 'i/'^iSzi). A finite J' breaks the degeneracy between the 
states, but as long as J' <^ J, the occupation of states with Si < 3/2 is energetically 
suppressed. These states can only be populated "virtually," leading to an effective 
Heisenberg interaction between the various effective S" = 3/2 states at each plaquette [8], 
i.e., 

Heff = G^Si- Sy (2) 

{id) 

— * 

Here, Si = {S^i, Syi, S^i) are spin 3/2 operators acting on the pseudospin states 
IS* = 3/2, Szi), and we have set h = 1. The interaction coefficient can be written 
as G = gJ'^/J, where (7 > is an antiferromagnetic-coupling constant that slowly 
varies as a function of J/U. Equation ^ explicitly shows the fragility of Nagaoka 
ferromagnetism, since a weak couphng among the plaquettes leads to a many-body 
ground state with antiferromagnetic correlations. 

To overcome this limitation, we consider a different initial configuration. Starting 
with four atoms per plaquette in the lowest orbital, we excite one of the atoms to a 
nondegenerate vibrational level [see Fig. |3]. This system is described by a two-band 
Hubbard Hamiltonian of the form 

{r,r'},a,n rn r 

— J ex X^ '^lra'^lra''^2rcT''^2ra (3) 
raa' 

which is characterized by on-site interactions between particles in the ground (f/n = U) 
and excited {U22 = Ue) bands, tunneling in the lower (Ji = J) and upper (J2 = Je) 
bands, and direct (V) and exchange (Jex) interactions between the two bands. In the 
present implementation, V = J ex- In Eq. (|3]), we have neglected terms that transfer 
atoms between bands, since they are energetically suppressed. The energy splitting 
between them has been omitted in our rotating frame. 

A single plaquette with three atoms in the lowest band and the fourth in the second 
band exhibits a crossing between an S=l and an S=2 state (ferromagnetic state) at a 
value of U that is smaller than Ut- Figure |4|^a) shows the low energy behavior of this 
system. For the parameters used in Fig.|4[a), the crossing occurs at f/ ~ 6 J. In general, 
U depends on J, Je, f/, and Jex- Figure |4[b) analyzes the existence of a ferromagnetic 
ground state as a function ofU/J and aj = Je/ J. For this analysis, we assume that the 
interaction terms are proportional, i.e., Jex = oij^^U . Different shaded regions display 
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Figure 3. Schematic representation of two coupled double-band plaquettes. Solid 
circles represent occupied orbitals. 



the ferromagnetic regions for different aj^^, and the dashed curves their corresponding 
critical values U . As both aj and a^x increase, the U values decrease, extending the 
ferromagnetic region. This behavior is consistent with our physical picture that both 
double exchange processes and the Hund's rule coupling stabilize ferromagnetism. 
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Figure 4. (a) Energies of a plaquette as a function oiU / J and the scattering length 
in Bohr radii Uq. The parameters that characterize the Hamiltonian [Eq. ([3|] are 
obtained for a superlattice constructed with a short-wavelength laser of = 765 nm 
that characterizes the short-lattice recoil energy Er — h'^/{2mX'^). The energies of the 
figure corresponds to Vi = 20 Er and Vs = 7.5 Er- Inset: Zoom in of the spectrum at 
the ferromagnetic crossing. Solid lines are the spectrum for zero gradient field while 
dashed lines are the spectrum at finite gradient field, (b) Ferromagnetism in a two- 
band plaquette. Shaded areas correspond to a ferromagnetic ground state. Dashed 
curves correspond (from top to bottom) to aj^^ = 0.2, 0.4, 0.6, and 0.8. 



The mobile atoms in the excited band are also expected to stabilize, via double- 
exchange processes, the ferromagnetic phase when a weak tunneling between plaquettes 
(J' and Jg) is allowed. The stabilization mechanism relies on the preservation of the 
spin when an atom hops, and on the energy penalty of 2Jex when ground and excited 
atoms form a singlet instead of a triplet at a given site. Only when the spins of 
adjacent plaquettes are fully aligned, the mobile atoms are free to hop. We confirm the 
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stabilization of the ferromagnetic correlations by studying the weakly coupled regime, 
which can be described again by an effective Heisenberg Hamiltonian as in Eq. Q, but 
now between the S = 2 states at each plaquette. 




Figure 5. (a) Asymptotic coefficients g°° and g'^ as functions of aj. The solid 
curve corresponds to g, and the dotted curve corresponds to g^. The dashed curve 
corresponds to g^ when ajj^ = 0. (b) Phase diagram for J' = 0. The shaded regions 
correspond to the ferromagnetic regime for different aj values. From left to right, 
aj = 0.5, aj = 1.1, aj = 2. 



The coupling coefficient G can be obtained by considering virtual processes in 
which one atom from the ground (excited) band hops from one plaquette to the other 
and returns to the original configuration. In practice, we extract G from the analysis 
of the low-energy spectrum of two weakly coupled plaquettes [see Fig. |3] . G depends 
in a nontrivial way on the interaction and kinetic energy parameters and it is given 
by G = gJ'^/J + geJ'e/Je for J', Jg ^ J,Je- Here g and ge are proportional to 
G when the interplaquette tunneling is allowed only in the ground or excited band, 
respectively. Both depend on J, Jg, U , and J ex] and g^ also depends on U^. In the 
effective model, the existence of ferromagnetism can be directly deduced from the sign 
of G. To further understand the robustness of the ferromagnetic phase, we consider 
the general case for which all the interaction terms are proportional to each other, 
Ue = aueU, V = Jex = ctJe^U, and the kinetic terms are related as Je = ajJ. For large 
interaction values {U — )■ oo), the coefficients g and ge approach to their asymptotic 
values g°° and g'^ that only depend on aj (independent of a j^^ and au^ so long as 
aj^^ > and au^ > 0). Figure [s] (a) shows that both g°° and g"^ become negative 
for a large parameter regime, confirming the robustness of the ferromagnetic phase. If 
Ue = [dashed curve in Fig. [s] (a)], g"^ becomes positive for an important region of aj 
values. This finding is consistent with Ref. [H] where it is pointed out that a nonzero 
interaction between atoms in the excited band {Ue > 0) can be crucial for the transition 
to a ferromagnetic ground state. 

Figure [s] (b) shows the phase diagram as a function of au^ and U/ J for three 
different aj = 0.5, 1.1, 2 values, with the assumption that J' = 0. For this study, we 
set aj^^ = 0.4, which is a typical value for optical superlattices. In general, we observe 
that as aj^^ increases, ferromagnetism becomes more favorable. The transition to a 
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ferromagnetic ground state depends strongly on a j and, for the J' = case, low values 
of a J enhance ferromagnetism. For Jg = 0, the phase diagram is independent of au^, 
and higher mobility of the excited atoms, large aj, favors ferromagnetic correlations. 
For example, for aj = 2 the critical value Uc is ~ 147 J. It increases to Uc ~ 488J for 
a J = 1.1, and for aj = 0.5, no ferromagnetic phase is observed. 




Figure 6. Lowest energies as a function of U/J and tiie scattering lengtli of two 
weakly coupled plaquettes (V/ = 4Er and Vs — 5.5Er) with six particles in the lowest 
band and two in the excited band and a total Sz — 0. Circles correspond to exact 
numerical calculations, and lines correspond to the effective Hamiltonian [Eq. ([2|] 
description. Top grid: scattering length values in Bohr radii for ^Li parameters. Inset: 
The G coefficient as a function of U/J for a plaquette with Vi = AEr and Vs = 5.5Er 
(dashed curve). For this case, the tunneling is J sa 0.085£^r- The transition occurs at 
U/J w 100. 



However, in realistic experimental setups it is not possible to explore the complete 
parameter space. For standard superlattice geometries, in which only the ratio Vi/Vg 
controls the relation between the different tunneling parameters, we find that a favorable 
ferromagnetic scenario takes place in the regime when J = Jg and J' = Jg. We 
can achieve this scenario by loading the atoms in the lowest vibrational state of a 
2D plaquette array (tight confinement along z) creating a 2D lattice in the x — y 
direction and and then exciting one atom per plaquette to the first vibrational 
orbital along the z direction. The atoms are initially all in the ground-state orbital 
$g(r) = 0o(a;)(/>o(2/)0o(2^), and one of the atoms in each plaquette is excited to the 
$e(r) = 0o(a^)0o(y)0i(-2) vibrational state. The tunneling along the x and y directions 
is independent of the z-orbital {(j)o{z) or (j)i{z)) and, therefore, J = Je and J' = Jg. 
The tunneling in the z direction depends on the z-orbital but it is negligible in the 2D 
geometry in consideration. 

Figure [6] presents this case for realistic ^Li experimental parameters (aj^^ = 0.5 and 
au^ = 0.7) where we observe a change in the sign of the coupling coefficient G (see inset) 
from positive to negative at f// J ~ 100, a value of interaction accessible with a Feshbach 
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resonance. Even though for the perturbative regime under consideration, where the 
effective Hamiltonian is vahd, G is small ( of the order of few Hz), it is measurable with 
state-of-the-art technology, as demonstrated in recent experiments where superexchange 
interactions as low as 5 Hz have been resolved [T3] . 

For preparing and probing ferromagnetism in the isolated two-band plaquette array, 
we propose to start by loading atoms into a band insulator in a deep 3D lattice 116] . 
By applying a double-well superlattice along the z-direction and manipulating the 
double-well bias, one can excite one of the four atoms in each double well to the first 
vibrational state via interaction blockade [17]. Then, by slowly merging the double wells 
along the 2;-direction and splitting them along the x — y directions, one can load the 
desired plaquette configuration with three atoms in the lowest and the fourth in the first- 
excited vibrational state along z. Alternatively, spatially selective two-photon Raman 
pulses can be applied to excite one atom in each plaquette to the desired vibrational 
state [m [191 ISO]- This procedure will lead to plaquettes with 5" = 0. The state 
5* = can be adiabatically converted into 5* = 1 by applying a nonlinear magnetic field 
gradientjH along z and tuning the interactions from f/ < to f/ > 0. The transition 
between S = 1 and 5* = 2 states in a single two-band plaquette can be probed using the 
same experimental techniques proposed for the single band plaquette. A magnetic field 
gradient couples the S = 1 and S = 2 states [see inset of Fig. |4|^a)] and this coupling 
can be used to probe the transition with spin-imbalance measurements or band mapping 
techniques. 

Interaction induced ferromagnetism in a few coupled plaquette array could also be 
observed by varying U/J in the presence of a magnetic field gradient. Since our initially 
prepared state is a band insulator with S = the magnetic field gradient has to be 
large enough to couple S = with S = Smax across the transition. After this procedure 
is applied the ferromagnetic nature of the ground state can be inferred in the applied 
magnetic-field gradient by measuring the local magnetization of the system [211 ESI ES] . 
A linear magnetic-field gradient produces a perturbation in the effective Hamiltonian of 
the form Hp = i6EpSzi/h, where 6Ep is the average energy shift between consecutive 
plaquettes. In the ferromagnetic phase, the formation of a domain wall is expected. The 
domain-wall width will be determined by the dimensionless parameter zGS/ 6Ep, where z 
is the number of nearest-neighbor plaquettes, and S = 2. The measurement of this width 
can be used to extract G in the ferromagnetic regime. In the antiferromagnetic phase on 
the contrary, no domain wall will be formed, and the local Neel order parameter should 
vary smoothly. The onset of ferromagnetic correlations as the system is driven through 
the critical point should be signaled by a suppression in inelastic collisions, a minimum 
in kinetic energy, and a maximum in the size of the cloud. The latter signatures have 
been demonstrated to be useful smoking guns of a ferromagnetic transition in recent 
experiments carried on in fermionic gases without a lattice (See Ref. [21] and references 
therein). 

§ Our calculations show that a simple linear gradient does not couple the 5 = 1 and S — states. 



Probing Interaction- Induced Ferromagnetism in Optical Superlattices 



11 



4. Conclusions 

In summary, we have proposed a controllable and experimentally realizable scheme to 
study interaction-induced ferromagnetism in ultracold atoms. Our method exploits the 
advantage offered by these systems to divide the full lattice into plaquettes. 

We showed that a plaquette loaded with three fermionic atoms is a promising 
set-up to experimentally observe Nagaoka ferromagnetism for the first time. We also 
used the plaquettes as the fundamental building blocks to create interaction induced 
ferromagnetism at longer length scales. We analyzed the system of weakly coupled 
two-band plaquettes and demonstrated that in the limit where the interplaquette 
coupling can be treated perturbatively, the system maps out into an effective Heisenberg 
Hamiltonian with a coupling constant G which changes sign from positive to negative as 
interactions are increased. The change in sign is a manifestation of an antiferromagnetic 
to ferromagnetic transition. 

In the perturbative regime, where our analysis is valid, the coupling parameter 
G is small. Consequently, it will be experimentally challenging to adiabatically reach 
the ferromagnetic ground state for large number of coupled plaquettes by increasing 
interactions. In this situation, the required temperature would be smaller than 
G ~ 10~^J. Nevertheless, exact diagonalization in a two-plaquette array confirmed the 
persistence of itinerant ferromagnetism beyond the weakly coupling regime. Actually, 
we found excellent agreement between our perturbative Hamiltonian and the many- 
body spectrum even at values of J' / J as high as 1/4, as shown in the low energy 
spectrum presented in Fig. [6j Our two-plaquette results seem to be consistent with 
variational Monte Carlo [21] and dynamical mean-field [2S] predictions that have found 
ferromagnetic phases in the two-band generic square-lattice Hubbard model. The 
actual stability of the ferromagnetic phase in larger plaquette arrays and stronger 
interplaquette couplings will need however to be resolved ultimately by experiments. 
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